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ABSTRACT 



The pressure-driven inertial collapse of a cylindrical void in an inviscid liquid is an 
integrable, Hamiltonian system that forms a finite-time singularity as the radius of 
the void collapses to zero. Here it is shown that when the natural cylindrical sym- 
metry of the void is perturbed azimuthally, the perturbation modes neither grow 
nor decay, but instead cause constant amplitude vibrations about the leading-order 
symmetric collapse. Though the amplitudes are frozen in time, they grow relative to 
the mean radius which is collapsing to zero, eventually overtaking the leading-order 
symmetric implosion. Including weak viscous dissipation destroys the integr ability of 
the underlying symmetric implosion, and the effect on the stability spectrum is that 
short-wavelength disturbances are now erased as the implosion proceeds. Introducing 
a weak rotational flow component to the symmetric implosion dynamics causes the vi- 
brating shapes to spin as the mean radius collapses. The above theoretical scenario is 
compared to a closely related experimental realization of void implosion: the discon- 
nection of an air bubble from an underwater nozzle. There, the thin neck connecting 
the bubble to the nozzle implodes primarily radially inward and disconnects. Recent 
experiments were able to induce vibrations of the neck shape by releasing the bubble 
from a slot-shaped nozzle. The frequency and amplitude of the observed vibrations 
are consistent with the theoretical prediction once surface tension effects are taken 
into account. 
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1.1 A. Air bubble released from underwater nozzle, images from N. Keim 
(U. Chicago) B. Close-up images of radial implosion of neck as it dis- 
connects. C. Nozzle is tilted by two degrees and shapes before and 
after disconnection show signs of the asymmetry. Two tiny side-by- 
side satellite bubbles are left behind, evidence of the inherent memory 
in this system. Nozzle is tilted towards viewer in the 'front' view, and 
to the left side in the 'side' view 4 

3.1 Numerical solutions to (3.4) for vibrational mode amplitudes in inertial 
hmit, for modes n — 2 (dashed line), n — b (dotted line). The ampli- 
tudes grow weakly as the mean radius decreases, like (ln(i?oo/i?))''"/^, 
predicted from solution to (3.5) using (thick solid line). Mode ampli- 
tudes are rcscaled by their initial size feg ^"^^ the mean radius R is 
rescaled by its initial size Rq, and Rqq/Rq = 100 18 

3.2 Illustration of void shapes for superposition of vibrational modes n = 2 
and 3, with dynamics given by Eq. 3.7 in text. The shapes arc rcscaled 
by the mean radius which is collapsing to zero (left to right). As the 
mean radius collapses, the shapes become more distorted. Because the 

n = 2 and 3 modes are not vibrating in phase, complex shapes ensue. 20 

3.3 Example of disturbance to the pressure and velocity fields outside the 
surface due to an n = 4 shape distortion (one quadrant shown). For 
the pressure contours, inside the void the disturbance pressure is zero 
(medium grey) and outside the pressure has regions of positive (lighter) 
and negative (darker) pressure. The surface location is highlighted by 
the thick black line in the lower plot for the velocity. At this instant 
in time, the mean radius R — lO/im and has an implosion speed R — 
—0.7fira/ fis. For clarity, a large perturbation amplitude is shown with 

64 = 8/xm and 64 = — 0.1/xm//xs 22 

4.1 Amplitude decay due to viscous dissipation at the interface. Here, 
Rqq/Rq = 10'^ and the dynamic Reynolds number for the flow is 
Q/u = 10"^ when R = Rq. The amplitudes are rescaled by their initial 
values 60 and the mean radius is rescaled by its initial value Rq. (a) 
Representative vibration for n = 2 (thin line) shows that viscous dissi- 
pation causes a decrease in the amplitude (envelope is thick line) and 
frequency of vibration, (b) Envelope of vibrational mode amplitudes 
for n = 2 — 10 (upper- to lower-most). Short wavelength (high n) 
disturbances are efficiently smoothed out by the presence of viscosity, 
(c) Amplitude spectrum taken at R/Rq — 10~^ (open circles) and at 
R/Rq = 10-5 (solid circles) 27 
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4.2 Amplitude decay due to viscous dissipation at the interface. Here, 
Rqo/Rq = 10^ and the dynamic Reynolds number for the flow is 
Q/v — 10^ when R — Rq. The amphtudes are rescaled by their initial 
values foo ^iid the mean radius is rescaled by its initial value Rq. (a) 
Representative vibration for n = 2 (thin line) shows that viscous dissi- 
pation causes a decrease in the amplitude (envelope is thick line) and 
frequency of vibration, (b) Envelope of vibrational mode amplitudes 
for n = 2 — 10 (upper- to lower-most). Short wavelength (high n) 
disturbances are efficiently smoothed out by the presence of viscosity, 
(c) Amplitude spectrum taken at R/Rq — 10~^ (open circles) and at 
R/Rq = 10-5 (sohd circles) 28 

4.3 Coupled mode amplitudes with presence of swirl flow for n = 10, biQ 
(thick line) and ciq (thin line, displaced downwards for clarity). The 
initial conditions at R — Rq are biQ — ciq — Bq with zero derivative. 
The relative swirl flow strength is ug/ur ~ 10~^ (over the range of the 
plot it increases only to 3 ■ 10"^ as decreases). Roo/Rq = 10^ and 

the implosion energy is about 32 

4.4 Effect of surface tension on early (large R) vibrational mode behavior. 
Here, Rqc/Rq = 10 and E/tt^Rq — 10. The amplitudes are rescaled 
by their initial values bQ and the mean radius is rescaled by its initial 
value Rq. (a) Representative vibration for n = 10 (thin line) shows 
that surface tension causes an initial increase in the amplitude (enve- 
lope is thick line) and frequency of vibration. After the mean radius 
has decreased beyond R^^ ~ -Rq/IO^, the vibration settles in to the 
asymptotic inertial state with amplitudes and periods as in Fig. 3.1. 
(b) Envelope of vibrational mode amplitudes for n = 2 — 30 (lower- to 
upper-most), (c) Amplitude spectrum taken at R/Rq = 10"^ (open 
circles) and at R/Rq = 10~^ (solid circles) 35 

5.1 (a) Collapse of minimum radius for air bubble released from circular, 
4mm diameter nozzle (red circles) . Theoretical inertial implosion from 
(2.7) with an implosion energy per unit length of = 12 dyne (black 
line). Dashed hne shows slope of 1/2. (b) Experimental neck profiles 
of the same air bubble (blue diamonds) . Theoretical profiles calculated 
as described in text using conservation of energy for each isolated layer 
(lines). Profiles are at 4.5, 49, 108, 166, 270, and 372 /is before dis- 
connection (inner to outer curves). Experimental data provided by 
Nathan Keim (U. of Chicago) 39 



Experimentally observed vibrations of n = 2 mode by release from 1.6 
mm by 9.6 mm slot nozzle (red). Theoretical vibration calculated from 
(4.8) as described in text (solid black line). Inset: Experimental mean 
radius collapse (red) corresponds to an inertial implosion with E — 12 
dyne, and i?oo = 5mm (solid black). A power-law of 1/2 is shown by 
the dashed line for comparison 



CHAPTER 1 
INTRODUCTION 



One of the most fundamental questions in fluid mechanics, is how does a volume of 
fluid break into two? It happens all the time, all around us - for example a stream of 
water faUing from a faucet breaks into droplets, then these droplets splash as they hit 
the sink, sending off even more tiny satellite droplets. Despite the ubiquitous nature 
of the phenomena, we still do not have a complete understanding of the process for 
a generic combination of fluids [1, 2]. The question is even more interesting because 
when a fluid undergoes a topological transition, a singularity occurs in the governing 
Navier-Stokes equations, where physical quantities such as the liquid pressure or 
velocity diverge in a finite time [3, 4, 5]. The mathematical appeal of singularity 
formation and technological relevance of the generation of droplets and bubbles has 
spurred much interest in topological changes of fluid surfaces in past years. Many 
of the theoretical works have concentrated on characterizing the universal features of 
the approach to the singularity [6, 7, 8]. 

However, recent results from multiple experimental groups have indicated that 
the pressure-driven implosion of a gas cavity in an inviscid liquid is not universal and 
that small asymmetries present in the initial state are preserved. The retention of 
such disturbances has the potential to dramatically alter the final state at breakup. In 
these experiments, a cavity is created and then disconnected either by quasi-statically 
releasing a bubble from a nozzle submerged in a hquid (see Fig. 1.1) [9, 10, 11, 12, 13], 
by impacting a solid disc into a liquid bath [14, 15, 16], or by placing a bubble in a 
co-flowing stream [17]. The cavity forms following the bubble as it rises, following 
the disc as it passes into the bath, or in the middle of the bubble as it is pulled 
apart by the divergent flow. For all cases, the final dynamics before disconnection is 
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dominated by a rapid inertia-driven implosion, not by surface tension which drives 
the breakup of a hquid drop in air [6, 18, 19]. In the only detailed investigation of 
the effect of initial azimuthal asymmetry on the bubble pinch-off dynamics Keim et 
al [10] found that tilting the nozzle by 2° results in side-by-side satellite bubble 
formation and non-circular neck cross sections. Tilts of even less than a tenth of a 
degree result in visibly asymmetric neck profiles, demonstrating that the generic form 
of pinch-off realized in experiments is asymmetric. Despite the fact that azimuthal 
asymmetry is usually observed in the experiments, it has only recently begun to be 
explored theoretically [13]. 

The inertial implosion of a gas cavity belongs to a different class of dynamics which 
are not universal, but instead have a precise memory of the details of the initial state. 
To understand the origin of the observed memory, consider a theoretically simpler 
but closely related system - the implosion of a cylindrical void in an inviscid liq- 
uid. Integrable systems such as this have a perfect memory (with as many conserved 
quantities as independent variables), and when they are perturbed, the Kolmogorov- 
Arnold-Moser theorem predicts that the dynamics will have trajectories which, for 
the most part, closely follow those of the original integrable dynamics [20]. In this 
system, which is both integrable and singular, the underlying memory manifests itself 
in a particularly simple way - as the amplitudes of shape disturbances become frozen 
in time as the cavity implodes [13]. The constant-amplitude disturbances cause vi- 
brations of the void's cross section shape as the implosion proceeds, with frequencies 
proportional to 1/r, where r is the time until the radius closes to zero. Since the 
frequency diverges as r — > 0, the shapes chirp, vibrating faster and faster as the void 
closes. This presents a natural explanation why the experiments observe a variety of 
disconnection scenarios depending on the specific experimental conditions. Since the 
experiments have a dynamics close to an integrable dynamics, small perturbations 
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are retained as the hole closes. 

The first purpose of this paper is to investigate additional physical effects which 
are mathematically interesting and experimentally relevant. Including weak viscous 
dissipation in the model destroys the integr ability of the leading order implosion dy- 
namics, and including a non-zero circulation to the liquid fiow prevents the formation 
of a singularity. Both effects have consequences on how a disturbance on the void sur- 
face evolves over time, but within the experimentally accessible regime and for large 
wavelength (small mode number) disturbances do not qualitatively change the freez- 
ing and chirping behavior. The second purpose is to directly connect the theoretical 
model to the experimental realization for both the symmetric and perturbed cases. 
This requires including surface tension, which, though it becomes asymptotically ir- 
relevant compared to inertial forces, has lingering effects on the vibration frequency 
within the experimental regime. 

We begin by reviewing the cylindrically symmetric implosion dynamics and demon- 
strate that the dynamics can be recast into Hamiltonian form. We examine the 
model's implications for quasi-2D (nearly cylindrical) void collapse and show that 
though the integrable dynamics is fundamentally non-universal, it appears self-similar 
in this limit. Then we consider the stability of the cylindrically symmetric implosion 
when subjected to azimuthal shape perturbations and derive the mode amplitude 
freezing and chirping behavior, including asymptotic logarithmic corrections. Next, 
we investigate changes to the leading-order symmetric dynamics and vibrational mode 
behavior after taking into account viscous dissipation and introducing a rotational 
fiow. Finally, we discuss the relation between experiments for the cylindrically sym- 
metric and azimuthally perturbed cases, ending with a direct comparison between the 
theoretical prediction for a single vibrational mode and an experimentally induced 
vibration in the neck of a bubble released from an oblong nozzle. 




Figure 1.1: A. Air bubble released from underwater nozzle, images from N. Keim (U. 
Chicago) B. Close-up images of radial implosion of neck as it disconnects. C. Nozzle 
is tilted by two degrees and shapes before and after disconnection show signs of the 
asymmetry. Two tiny side-by-side satellite bubbles are left behind, evidence of the 
inherent memory in this system. Nozzle is tilted towards viewer in the 'front' view, 
and to the left side in the 'side' view. 



CHAPTER 2 
AXISYMMETRIC CAVITY IMPLOSION 

The implosion of an arbitrarily shaped cavity is a very complicated, fully 3D prob- 
lem. However, with additional geometric constraints the problem is simplified. The 
dynamics of a spherical cavity have been studied for over a century [21, 22, 23], and 
is an integral component in the understanding cavitation [24] and the phenomena of 
sonoluminescence [25]. The inertial collapse of a cylindrical cavity was first studied in 
relation to bubble production from underwater nozzles [26, 27] and has more recently 
been applied to cavity collapse in liquids after solid body impact [14] and cavities 
in granular media [28]. The inner gas dynamics can also play a role in the asymp- 
totic dynamics [29] . Here, we briefly reexamine the cylindrically symmetric implosion 
and its implications so that we can later characterize its stability when perturbed 
azimuthally. 

The problem under consideration is the inertial implosion of an infinite cylindrical 
vacuum cavity in an inviscid liquid. After the void is created it begins to implode 
inwards because the pressure inside the void is less than the pressure of the liquid 
outside. The inviscid fluid flow is governed by Euler equations: du/ dt + u • Vu — 
—VP/p), where u is the local velocity field, P the pressure, and p the liquid density. 
When incompressible, the continuity equation for volume conservation V ■ u = 
apphes. If the process starts from rest, the flow must also be irrotational {V xu — 0) 

— * 

for all times and thus we can also write the velocity asu — V$ where $ is the velocity 
potential which satisfies Laplace's equation V $ = 0, following from the continuity 
equation. Solving for the flow is just a matter of solving Laplace's equation, which is 
linear and time-independent. The difficulty in this problem arises via the non-linear 
dynamic boundary conditions (the normal stress balance on the interface) which 
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determines how the interface moves over time. 

To describe the implosion of the cyhndrical cavity, we naturally use circular cylin- 
drical coordinates, with the z-eods aligned along the axis of the cylinder, and the 
surface at r = R{t). The implosion dynamics is completely 2;- independent, and two- 
dimensional because we can consider the implosion of a circle of radius R{t) in a sheet 
of water. If no rotational flow uq is present, the solution from the Euler equations for 
the flow field is purely a convergent radial flow of the form Ur{r,t) = —Q{t)/r and 
velocity potential $ = —Q{t) ln(r). This corresponds to the flow of a point-sink with 
strength Q(t), which will be determined by the boundary conditions. The strength 
Q{t) is related to the areal flux of incoming water by dA/dt = 27tRR = 2T:Q{t). The 
first boundary condition on the motion of the interface is the kinematic condition 
that the interface moves with a velocity equal to that of the normal velocity of the 
liquid at the interface: 

Mr,t)\nit)--^-R, (2-1) 

showing that the flux Q{t) — —RR and the velocity fleld at any point in time are 
completely determined by the radius of the hole and it's velocity, Ur{r, t) — RR/r. 

To flnd the evolution of the radius R{t), Euler's equations can be integrated over 
the region of the liquid (from r = R{t) to Rqo) where Rqo is a far away distance 
where the radial flow has decayed or where the inner flow transitions to an outer 
flow of a different form. In 2D we must impose this length scale to keep the energy 
flnite because the potential does not decay as r — > 00. Experimentally, this scale 
will depend on the system size and set-up. Integrating Euler's equations yields the 
Bernoulli integral [23] for the pressure in the liquid P{r,t). It essentially constitutes 
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a conservation of energy for each fluid volume. In this case it yields 




(2.2) 



where Poo is the pressure at Rqo taken as constant. The normal stress balance at 
the imploding surface is given by Laplace's formula for the pressure jump across the 
interface P{R, t) — Pq — ^k, where Pq is the cavity pressure taken as constant, 7 is the 
line tension (analog of surface tension in 2D), and k is the curvature which is 1/R{t) 
for a circle. By evaluating the pressure in (2.2) at the surface, we can eliminate it via 
the normal stress balance and derive an evolution equation for the radius R{t): 



In (2.3) we have incorporated the fact that the pressure Pqo is higher than Pq by 
an amount AP. Ignoring the terms in (2.3) which are 0(P/Poo) or smaller and 
inserting the known form for $(t) yields the well known 2-D analog of the Rayleigh- 
Plesset equation 



where the dot denotes a time derivative [24] . 

Equation (2.4) is the second-order non-linear ordinary differential equation which 
determines the evolution of the radius R{t) as it collapses to zero. Any one dimen- 
sional ordinary differential equation of the form {q -\- f{q)(p' -\- g{q) = 0), like (2.4) 
above, can be integrated and recast as a Hamiltonian [30]. We are not surprised by 
this because in fact any free boundary problem with inviscid, irrotational, and incom- 
pressible flow is Hamiltonian by nature [31]. However, whereas the Hamiltonian for 
an arbitrary flow may be quite complicated, here the Hamiltonian takes a particularly 




(2.3) 



AP + ^ = _ ^(^2 ^ ln(Poo/i?), 



(2.4) 



simple form and represents the total energy per unit length of the system. The fact 
that the dynamical equation (2.4) can be integrated has been remarked before [12, 14] 
but the implications of the integrability had not been fully exploited in this context. 
The integrability of the leading order dynamics combined with the approach to a 
singularity (where external length and time-scales are irrelevant) is the origin of the 
unusual fiat mode amplitude stability spectrum observed in this system [13] and other 
focusing phenomena like spherical void collapse [24] or cylindrical [32] and spherical 
shock wave implosion [33]. 

The total energy is the sum of the kinetic energy of the water in motion and the 
potential energy associated with opening up the hole in the liquid: 

p2 

H{R, Pr)^E= + APrf + 7 27rit:, (2.5) 

where R is the void radius, Pji = M[R)dR/dt is the radial momentum of the implo- 
sion, M{R) = 2p7ri?2 ln{Roo/R) is the mass of liquid in motion. The Hamiltonian is 
conserved and thus the total implosion energy is a constant. The kinetic term is 
simply the total kinetic energy of the sheet of water in motion: 

K.E. = pj^ 4 rdrde = -M(i?)i?2 = (2.6) 

The potential energy due to the pressure APnR^ comes from the work that must 
be done against the external pressure to form the hole and the surface energy ■y2TTR 
is the energy cost for adding surface length to the hole. One can easily check that 
this formulation is equivalent to the 2-D Rayleigh-Plesset equation (2.4) by applying 
Hamilton's equations of motion (i? = dH/dPji, Pji = —dH/dR). One technical 
advantage of using the equation in Hamiltonian, or integrated form, is that now we 
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no longer need to solve a 2nd order ODE for R{t), but can simply integrate the 
equation for dR/dt, which in fact has known solutions [12, 14]. A more fundamental 
imphcation is that it is now clear that this system belongs to the class of integrable 
dynamics as it has one degree of freedom and one conserved quantity, the energy. 

2.1 Asymptotic energy balance and closure dynamics 

For a void which starts from rest, initially the energy is due completely to the potential 
contributions from the overpressure in the liquid and surface tension. From (2.5) we 
can see that as the void closes, — > 0, this energy is transferred into kinetic energy. 
The asymptotic energy balance is 

E = -M{R)if = 7rpln{Roo/R){RRf. (2.7) 

A singularity occurs because the closure rate R must diverge to compensate for the 
decreasing area of the void in order to maintain constant energy. Another way to say 
this is that since the mass of liquid in motion M{R) decreases to at the moment 
of disconnection, the closure rate R must diverge. Both surface tension and the 
pressure effects thus become negligible relative to inertia as the implosion proceeds 
and R diverges. 

The leading order dynamics in the inertially dominated regime can be found by 
examining the right hand side of (2.4) or equivalcntly the energy balance in Eq. (2.7). 
From the energy balance, we see that the areal flux Q — —RR — -sj E / ph\{Roo / R) 
is slowly varying as i? — > 0. The flux slowly approaches zero as less and less liquid 
is required to move in to fill the hole. Approximating the flux Q as constant, yields 
the power law scaling R = a{ti) — t)^/'^ flrst proposed by Longuet-Higgins et al in 
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a study of the sound emission from bubbles released underwater [26]. Here, tjj is 
the closure time, and a is related to the flux Q hy — 2Q with dimensions of 
diffusivity (lengths/time). The prefactor a depends on the experimental specifics. 
For a bubble released from an underwater nozzle the appropriate physical scaling is 
~ -yy^Rj^l/p, with R]\f the nozzle radius [10]. 

In (2.4) if we include the logarithmic term but treat it as nearly constant, we find 
an acceptable power law solution R{t) = a{ti, — t)^ , with a weakly varying apparent 
exponent, 

41n(i?oo/-R) 

The exponent (3 tends towards 1/2 very slowly as i? — )■ 0, but remains slightly greater 
than 1/2, consistent with experimental observations of exponents in the range 0.54 — 
0.57 [9, 14, 10, 11]. This is a refinement of the first estimate of power law of ^ = 1/2 
from [26] to include the weak logarithmic dependence on R, consistent with previous 
analyses [14, 29]. 

2.2 Axial variation 

The above analysis is for the radial implosion of a perfectly cylindrical void, or equiv- 
alently, the implosion of a circular void in a 2D sheet of inviscid liquid. To extend 
the results to make predictions about how the full profile of 3D cavity (hke the nearly 
cylindrical neck connecting an air bubble to an underwater nozzle) collapses we make 
a long-wavelength approximation for the neck radius so that the variation in the ra- 
dius satisfies dR/dz <^ 1. This approximation was first used in [26] with the leading 
order dynamics treated as quasi-2D, and follows the above analysis for a cyhnder but 
with z simply treated as a parameter. The flow at each height converges radially 
with flux Q, and is isolated from any coupling in the ^-direction (there is no energy 
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or momentum transfer along z). For simplicity, we first examine the result ignoring 
logarithmic corrections so that the flux is constant in time (corresponding to the 
R ~ ar^/"^ power law) and in z. 

We initialize the surface at some t = within the inertially dominated implosion 
regime. In the region near the neck minimum the shape is quadratic 

R[z, t = 0) = H{z) = Hq + cqz^. (2.9) 

To stay within the quasi-2D limit the aspect ratio is constrained to be small, Hq/L — 
^/cqHq <^ 1, where we have used the axial length scale L — -s/Hq/cq. Taking the 
approximation that the flux is constant, we can write the evolution of the profile as 

R{z,t)^a{mz) -t)^/'^, (2.10) 



with the only difference between the evolution of the radius at difference heights being 
the implosion time, ti,{z), which is set by the initial radius profile. The implosion time 
is the time when the radius at that height will collapse to zero as determined by the 
initial size (2.9): 

.w^(^)^. 

With nothing to break the up/down symmetry about 2; = 0, the height which starts 
with the minimum radius Hq, will be the height of the disconnection event &X,t — tif — 
^5(0). Writing the shape evolution in terms of the time until disconnection t — t^—t 
yields 

By proper rescaling of the axial and radial coordinates, the evolution has a self-similar 
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form. Taking ( = R/hj^in V = ^/hmin with /imin = -^(^ = ^1^)1 the shape in the 
rescaled coordinates is approximated by: 

C = (2coi^or;2 + i)l/2_ ^2.13) 

In the hmit h^i^ — >■ 0, the shapes becomes conical {R{z) — \/2cqHq z) as seen 
in the experiments (Fig. 1.1) and shown in [26]. The relevant axial length scale 
is proportional to h^i^, so that it decreases in pace with the radial scale and the 
aspect ratio of the shapes are maintained. Of importance is the dependence of the 
self-similar form on the initial aspect ratio \/cqHq. The final slope of the cones is 
directly dependent on the initial aspect ratio of the quadratic neck shape. An initially 
slender shape will remain slender and the long-wavelength approximation analysis is 
self-consistent in this respect. 

Accounting for the logarithmic corrections to the power law exponent /3 in Eq. 2.8, 
we find a self-similar form 

C = (&^^ + 1)^, (2.14) 

with 77 = -z/^™,vf s-iid ^ — cqHV^ ^ I ^- The region where the shape is approximated 
by this form is bounded by 



77 < 



\ ^0(^ - l)Cn 



^0 (2.15) 



which grows as /i^m shrinks to 0. This form assumes that nearby layers evolve 
with the same a and a good approximation for slowly varying shapes because 
a and (3 themselves have only logarithmic dependencies on the local radius R{z,t). 
Accounting for the weak ^-variation of a and (3, and solving the governing equations 
self-consistently replacing Rqo with the axial length scale yields slight corrections to 
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the exponent [8]. 

Now, with the logarithmic corrections included, the axial length scale is propor- 
tional to h^J]^l^ which decreases almost in pace, but a bit slower, than the radial scale. 
For a power law measured over a finite range of time we can write /imin ^ r^l'^^^ so 
that the axial scale oc r-*^/^"*^, where b measures the difference in the exponent from 
1/2, b = l/41n(i?oo/^mm)- Other words, if the radius is decreasing with a power 
slightly greater than 1/2 (as observed in the experiments), then the quasi-2D model 
predicts the axial length scale will decrease with a power equally less than 1/2. This 
result is qualitatively consistent with experimental measurements of the axial radius 
of curvature at the minimum [14, 11] which find powers less than but still close to 
1/2. 

It is worthwhile to note that although the profiles have a self-similar form, it 
does not mean that the dynamics here is universal. Though often times universal 
dynamics exhibit self-similarity, in this case the reverse is not true. The underlying 
dynamics here is in fact nearly integrable, with the implosion energy at every height 
E{z) conserved. 

When extending the 2D-cylindrically symmetric model to capture the evolution of 
axisymmetric 3D shapes, we must also modify the Laplace pressure contribution due 
to surface tension at the interface which appears in the normal stress balance. Now, 
the surface is no longer a curve in a plane, but has an axial curvature component. 
Because the surface is curved outwards in the axial direction, it works in opposition 
to the azimuthal inwards curvature, effectively reducing the Laplace pressure jump 
across the interface. Directly from above, we know that the axial curvature at the 
minimum will increase at roughly the same rate as the radial curvature. This fact 
will later allow us to simply include the effect of the axial curvature by reducing the 
Laplace pressure by a fractional amount when comparing to the experiments over a 
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limited range. We can write the axial radius of curvature Rax = X^mim where x is 
a constant determined by the initial aspect ratio from (2.13). Thus, a more realistic 
value of the Laplace pressure jump 

= l{l/hmin - ^/Rax) = 7(1 - ^/x)/hmin (2-16) 

can be used in (2.3). 

The quasi-2D model outlined above provides a simple way to interpret and analyze 
the experimental results available today. It contains all the ingredients necessary to 
explain experimental observations of the time dependence of cavity profiles during 
implosion. We will see it is more than sufficient for the purposes of the present paper, 
not to predict the complete details of the asymptotic approach to the singularity, but 
to examine the effect of azimuthal asymmetry on the implosion dynamics within the 
experimental regime. 



CHAPTER 3 
AZIMUTHAL PERTURBATIONS 

With the axisymmetric implosion dynamics in the hmit of an infinite cyhnder now 
determined by equation (2.4) or conservation of the Hamiltonian (2.5), we can begin 
to investigate the stabihty of the implosion under a slight azimuthal shape distortion 
which destroys the symmetry. Because the underlying leading order dynamics is in- 
tegrable and possesses a natural memory mechanism via conservation of the energy, 
we will see the stability spectrum takes a special form. No single wavelength per- 
turbation will grow fastest to overtake the others, as typically occurs in interfacial 
instabilities. Instead, all the modes will grow at the same rate, logarithmically slow 
relative to the rapid collapse of the mean radius. This unique property is the result 
of the combination of the approach towards a singularity and the integrable dynam- 
ics [13]. Since all modes grow identically, the initial spectrum of mode amplitudes 
for a generic distortion is frozen in time as the singularity is approached, encoding 
an infinite amount of information about the initial state in form of vibrational modes 
about the radial implosion. 

Representing a distortion to the shape as the sum over the Fourier modes cos(n^) 
shifts the surface to he at 

r = S{e, t) = R{t) + J2 bn{t) cos(n^), (3.1) 

n 

where the bn{t) are the mode amplitudes and assumed to be small relative to mean 
radius R{t) {hn{i)/R{t) <^ 1). The velocity potential and velocity field in the liquid 
depend on the boundary conditions on the surface, so distorting the surface causes a 
disturbance fiow in the surrounding liquid. The leading order and disturbance fiow 
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have a total velocity potential of the form: 
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$(r, 0, t) = Q{t) ln(r) + J] dn{t)r-'' cos(n^). (3.2) 



n 



The kinematic condition at the surface (2.1) at 0{an/R) determines the coefficients 

dn{t) = — ibn + bn-\. (3.3) 

Note that here, the azimuthal components uq and hq caused by the disturbance 
did not enter into the kinematic equation (2.1) because they are both of 0{bn), 
thus making their product 0{b'^). Applying the BernouUi integral (2.3), and the 
Laplace formula for the pressure jump across the surface, also to 0{bn) yields the 
evolution equation for the mode amplitudes with coefficients which depend on the 
radial implosion dynamics: 



6n + I [|(l-n) |6n = 0. (3.4) 



We can take advantage of the fact that we have a conserved energy to use (2.4) and 
(2.5) to rewrite R and R in terms of E and R. This, along with using d/di = Rd/dR 
allows us to rewrite (3.4) so that the independent variable is R rather than t. Taking 
the limit that the implosion is inertially dominated and the radius has decreased far 
away from the system size Rqq, yields 

bn + (^ + + ^1^ = 0, (3.5) 

where ' denotes a derivative with respect to R. 

Immediately we can clearly see the amplitude freezing and vibrational behavior 
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by taking the limit ln(i?oo/-R) ^ in (3-5), which corresponds to a radial implosion 
follow it! = a.T^I'^ . Then, the governing equation takes the particularly simple equidi- 
mensional form 

&n + | + (^-l)^ = 0, (3.6) 



with solutions for bn{R) of the form bn oc i^^Vn-l^ These solutions neither grow 
nor decay because the eigenvalues are entirely imaginary, but instead oscillate with 
constant amplitude 6„ q n-dependent frequency. The real part of the vibration 

can be written 

bn{R) = bnfi cos{Vn - 1 In + 0„^o), (3-7) 

with (pnfi an arbitrary phase shift defined by the initial conditions. The frequency 
of the vibrations, y/n — 1/R is n-dependent, with higher modes vibrating faster. All 
the vibrations become faster and faster, or "chirp", as i? — )■ due to the divergent 
1/R dependence. This is a direct consequence of the singularity formation, where no 
other length or time-scales other than the closure radius R and time r are involved in 
the dynamics. Thus, from simple scaling the frequency is forced to be proportional 
to 1/R. 

To check that including the logarithmic corrections does not qualitatively change 
the vibration freezing and chirping characteristics, we analyze (3.5) by introducing a 
second slowly varying length scale defined as y — ln{RoQ/R). Performing a standard 
multiple scale analysis [34] we find that the modes are primarily oscillatory with am- 
plitudes that grow extremely weakly and uniformly, as (ln(i?oo/-R))^^^- This growth 
is shown in Fig. 3.1 and envelopes the numerical solutions to (3.4) for n = 2 and 5. 

The growth can also be shown by transforming (3.5) into the y variable: 

^_l^ + (,_l)fl_ 1^6^ = 0. (3.8) 
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10"''° 10"^ 10"^ 10""^ 10"^ 10° 

R/Rq 



Figure 3.1: Numerical solutions to (3.4) for vibrational mode amplitudes in inertial 
limit, for modes n — 2 (dashed line), n = 5 (dotted line). The amplitudes grow 
weakly as the mean radius decreases, like (ln(i?oo/-R))"'^^^, predicted from solution to 
(3.5) using (thick solid line). Mode amplitudes are rescaled by their initial size 6o 
and the mean radius R is rescaled by its initial size Rq, and Rqq/Rq — 100. 
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Solutions to (3.8) are of the form bn{y) — y3/2g-iv^n^ ^f{y) with f{y) satisfying 

2y——?r + —(5 — 4:iy/n — ly) + f(l-n- 5iV^^) = 0. (3.9) 
dy dy 

Again transforming to the complex variable z = 2i\/n — 1 y we arrive at 

which is of the form of the confluent hyper-geometric differential (Kummer's) equa- 
tion. The solution for j{y) when Re(z) 1 (or y 1, Rqq/R 1) decay like 
/ oc z~^l^. Taken together with the above form for hn{y) yields a weak growth in 
the amplitude like 6„(y) oc |/3/2-5/4 _ ^1/4 ^^(^^ ^ {\n{Roo / R))^ > ^ ■ 

The amplitude growth is so weak (in Fig. 3.1, increasing only by a factor of 1.5 
as the mean radius deceases by 10 orders of magnitude) that the amplitudes are 
effectively constant. More importantly, because all modes grow within the same 
n-independent envelope, an initial mode distribution will still be preserved as the 
implosion proceeds. Numerical solutions to (3.5) for n = 2 and 5 are shown in Fig. 
3.1. The weak growth follows the analytic prediction of ln(i?oo/-R)^^^- 

Considering the vibrational modes relative to the leading-order radial implosion 
the shape becomes more distorted as R decreases to 0. Though the amphtudes are 
constant, the relative amplitudes hn/R grow like r"-*^/^, causing the shape to become 
increasingly more distorted. Taking the linear stability results to the point when 
bn — 0{R) predicts the development of cusp-hke structures on the interface. The 
shapes, especially when two or more modes are present can become quite complex 
as the vibrations occur with n-dependent frequencies. An illustration of shapes that 
occur as the void collapses for a combination of modes n — 2 and 3 is shown in 
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Figure 3.2: Illustration of void shapes for superposition of vibrational modes n = 2 
and 3, with dynamics given by Eq. 3.7 in text. The shapes are rescaled by the mean 
radius which is collapsing to zero (left to right). As the mean radius collapses, the 
shapes become more distorted. Because the n — 2 and 3 modes are not vibrating in 
phase, complex shapes ensue. 

Fig. 3.2. 

At some point (once the mean radius R has collapsed to the size of the largest 
mode amplitude present) the approximation bn/R ^ 1 breaks down and non-linear 
effects cannot be ignored. Recent numerical results utilizing a conformal mapping 
method have explored the non-linear regime that follows and found a rich landscape 
of different types of pinch-off depending sensitively on the initial conditions [35]. 

These vibrational modes are purely inertial, unlike most vibrations or waves on 
fluid interfaces. For example, a jet of water issued from an elliptic shaped nozzle will 
also undergo vibrations of the cross section shape, but this effect is due to a competi- 
tion between inertia and surface tension [36] . The inertial shape vibrations here occur 
because once the flow and pressure flelds are perturbed, there is no mechanism for the 
extra energy associated to be dissipated in this closed and conservative system. The 
extra energy could go into or out of the main radial implosion, but in this situation 
it remains in the perturbation leaving the radial dynamics unchanged. 

To help visualize the dynamics of the inertial vibrations, we can track how the 
fluid moves in and out of the wells around the circle. The corresponding corrections 
to the pressure held and velocity fleld (after subtracting out the leading order radial 
implosion) can be calculated using the potential from (3.2) with the appropriate 
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coefficients from (3.3) given the evolution of the mode amphtude from (3.4). The 
disturbance pressure field P(r, 9) is found from evaluating the Bernoulli integral (2.3) 
at Oipn) and yields: 



For an n = 4 distortion the disturbance pressure and velocity fields are shown in 
Fig. 3. Physically, the pressure builds up in the wells, leading to a disturbance flow 
outwards towards the peaks. The pressure and velocity field oscillate, but are out of 
phase, so that as the shape rounds up to a circle, the pressure is uniform around the 
circle (in 9) but flow continues, and creates a well. The pressure in the well begins 
to build up again (lightest region in Fig. 3) , and causes a reversal of flow out of the 
wells. 

Thus far we have shown that when a small shape distortion is introduced, an 
imploding cylindrical void will undergo inertial shape vibrations about the radial 
implosion dynamics. Within the linear regime, these vibrations have constant ampli- 
tudes but grow relative to the shrinking mean radius. This scenario naturally provides 
an explanation for why experiments observe a wide variety of different disconnection 
events — the final moments at disconnection are sensitively linked to the earliest 
moments in the implosion process. Below we show that including additional physical 
effects into the model do not destroy this unique behavior within the experimental 
regime. 




(3.11) 



where Sn — hn-\- hnR/R- The disturbance velocity field is V$ at 0{bn)- 




(3.12) 
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Figure 3.3: Example of disturbance to the pressure and velocity fields outside the 
surface due to an n = 4 shape distortion (one quadrant shown). For the pressure 
contours, inside the void the disturbance pressure is zero (medium grey) and outside 
the pressure has regions of positive (lighter) and negative (darker) pressure. The 
surface location is highlighted by the thick black line in the lower plot for the velocity. 
At this instant in time, the mean radius R = 70/im and has an implosion speed 
R = — 0.7/im//is. For clarity, a large perturbation amplitude is shown with 64 = 8/im 
and 64 = — 0.1/im//is. 



CHAPTER 4 
CONSEQUENCES OF ADDITIONAL PHYSICAL 

EFFECTS 

The inertial shape vibrations were analyzed in a speciahzed hmit of flow and fluid 
properties. It was assumed that the leading order dynamics was a purely radial implo- 
sion, that the liquid was completely inviscid, and was analyzed in the far asymptotic 
limit when the dynamics was completely inertially dominated. For the purely radial 
implosion in the inviscid liquid, the dynamics was shown to be integrable and evolve 
towards a singularity. Below we will test the consequences on the vibrational modes 
of destroying these two properties, in ways that are physically relevant. 

This is done by first including viscous dissipation in the external liquid in the 
simplest manner possible, which destroys the integrability because the collapse energy 
is no longer conserved. Second, we remove the singularity by introducing a weak 
rotational fiow into the external liquid which prevents the hole from collapsing to a 
point and forming a finite time singularity. Both of these effects are physically relevant 
as well - no experiment uses a perfectly inviscid liquid or can completely eliminate the 
presence of a weak swirling fiow. We will find that within the experimental regime, 
the inertial shape vibration behavior described in the preceding section will be robust 
against these effects, for the small mode numbers (large wavelength disturbances) 
explored in the experiments. Finally, even though surface tension becomes irrelevant 
asymptotically, we account for capillary effects and find the vibration frequency is 
increased in the experimental regime before the vibrations settle in to their asymptotic 
behavior. 
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4.1 Viscous dissipation 

To include viscous dissipation we must first consider the effect on tlie leading-order 
radial implosion dynamics, and then consequently the effect on the vibrational mode 
behavior. When viscous dissipation is included in the model, the implosion energy 
E is no longer strictly conserved and the leading-order symmetric collapse dynamics 
are no longer integrable and Hamiltonian. However, Eq. (2.4) still holds once the 
viscous contribution to the normal stress balance is taken into account. For the 
leading-order radial flow, viscous dissipation does not alter the form of the equations 
governing fluid flow but it does contribute to the normal stress balance on the surface: 
p{R) = Po — + 2fidur/dr, where fi is the liquid viscosity. Since dur/dr — —R/R 
on the surface, the governing equation (2.4) is then modified to 

Ap + ^ + 2/i| = _ p^jf ^ ^jj) HRoo/R). (4.1) 

Comparing the inertial and viscous terms we see that viscous dissipation will be 
relevant to the process when RR/v — Q/v — 0(1) (with the kinematic viscosity 
v — n/p)i i.e., when the dynamic Reynolds number is 0{1). If the implosion be- 
gins in an inertially dominated regime where E ^constant, this occurs very late 
in the implosion, at a length scale beyond the reach of a continuum theory, seen 
from Q — E / t: p\n{RoQ / R) . For the air cavity in water experiments the dynamic 
Reynolds numbers are typically on the order of 10 — 10 , which means that within 
the experimental regime, the radial collapse is inertially dominated and E ^constant 
will produce the observed radial implosion dynamics. 

Once we consider the shape perturbations, the velocity field is no longer purely 
radial, but has an azimuthal {9) component. To fully include the effect of a small 
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viscosity in this problem would require a complete change in the solution method. 
The presence of viscosity means that vorticity can be generated and therefore the flow 
is no longer truly potential. However, if the vorticity is confined to a very shallow 
boundary layer near the surface where it is generated, we can again approximate the 
viscous effects as only contributing to the normal stress balance. Then, the stabihty 
analysis yields a revised equation for bji{t): 



From (4.2) we can see that viscous effects can now play an earlier role for high 
mode numbers, because the viscous terms have prefactors 0{n^) and will quickly 
begin to dominate. Figs. 4.1 and 4.2 plot results of numerical solutions to (4.2) for 
two values of dynamic Reynolds number, 10"^ and 10"^ (which over the range of the 
data change by only a factor of about 2). Figs. 4.1 and 4.2 (A) plot the vibration 
for a single mode n — 2, showing a decay in the amplitude relative to inviscid case 
(see Fig. 3.1), but the chirping behavior of the frequency is unchanged. In Figs. 4.1 
and 4.2 (B), the envelopes for modes n = 2 — 10 are plotted, and clearly show the 
n-dependence of the rate of decay. All mode amplitudes experience decay due to 
viscous dissipation, but the short wavelength (high n) feel the strongest effect. The 
reason for the apparent amplitude growth in Fig. 4.1 B for Q/v — 10^ is simply that 
the inertially dominated vibrations have weak growth (as in Fig. 3.1) which acts until 
viscous dissipation takes over as i? — ?> 0. For the smaller Q/u = 10"^ in Fig. 4.2, the 
weak growth is completely washed out by a rapid decay in the amplitude. 

To quantify the decay, we plot the spectrum of the mode amplitudes as a function 




(4.2) 
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of mode number n at a given mean radius. We choose two values of the mean radius 
at which to evaluate the spectrum, R/ Rq — 10~^ and 10~^ (so the mean radius 
collapses 2 or 5 decades from its initial value). The results are shown in Figs. 4.1 
and 4.2 (C). The choice of the mean radius at which to evaluate the spectrum has 
no qualitative effect on the results. Both spectrums show rapid decay for high mode 
numbers, with more efficient decay for the lower dynamic Reynolds number set. 

Overall, viscous dissipation does not effect the radial, symmetric implosion until 
exponentially small length-scales, but it has a dramatic effect on the vibrational mode 
behavior. The chirping characteristic of the vibrations persists, but the freezing of 
the mode amplitudes is lost and a rapid decay of the mode amplitudes occurs for high 
n and low dynamics Reynolds number. 

To connect this result back to the experiments, we note that in the bubble breakup 
experiments, surface tension is effective in smoothing out the surface and removing 
any higher mode distortions which may be present, because the early dynamics is 
completely dominated by surface tension before the implosion begins [9, 10, 11]. But, 
if a short wavelength disturbance was present, viscous dissipation would be very 
efficient at wiping it out (Fig. 4.2). So in these experiments we expect only low mode 
number distortions to be present. 

The solid body impact experiments do not have such a surface tension dominated 
regime and thus short wavelength disturbances may be present in the initial shape. 
With the faster radial implosion and higher dynamics Reynolds number in the exper- 
imental regime, viscous dissipation would cause decay of the highest modes but it is 
less effective here than in the bubble experiments so some higher modes could still be 
present, which could be part of the reason why the intricate textures and swarms of 
satellite bubbles are observed in those experiments [14]. 
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Figure 4.1: Amplitude decay due to viscous dissipation at the interface. Here, 
Roo/Rq = 10"^ and the dynamic Reynolds number for the flow is Q/u = 10"^ when 
R = Rq. The amplitudes are rescaled by their initial values 69 the mean radius 
is rescaled by its initial value Rq. (a) Representative vibration for n = 2 (thin line) 
shows that viscous dissipation causes a decrease in the amplitude (envelope is thick 
line) and frequency of vibration, (b) Envelope of vibrational mode amplitudes for 
n = 2 — 10 (upper- to lower-most). Short wavelength (high n) disturbances are effi- 
ciently smoothed out by the presence of viscosity, (c) Amplitude spectrum taken at 



R/Rq = 10 ^ (open circles) and at R/ Rq = 10 ^ (sohd circles) 
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Figure 4.2: Amplitude decay due to viscous dissipation at the interface. Here, 
Roo/Rq = 10"^ and the dynamic Reynolds number for the flow is Q/u = 10"^ when 
R = Rq. The amplitudes are rescaled by their initial values 69 the mean radius 
is rescaled by its initial value Rq. (a) Representative vibration for n = 2 (thin line) 
shows that viscous dissipation causes a decrease in the amplitude (envelope is thick 
line) and frequency of vibration, (b) Envelope of vibrational mode amplitudes for 
n = 2 — 10 (upper- to lower-most). Short wavelength (high n) disturbances are effi- 
ciently smoothed out by the presence of viscosity, (c) Amplitude spectrum taken at 
R/Rq — 10~^ (open circles) and at R/ Rq — 10~^ (sohd circles). 
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4.2 Rotational motion 
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Introducing a small rotational flow into the external liquid provides a clever way to 
remove the singularity formation that occurs as the radius of the void collapses to zero. 
This does not break the cylindrical symmetry, but causes the liquid to swirl around 
as the hole collapses. A consequence of the rotation is that finite amount of energy 
is always stored in the rotation flow because of angular momentum conservation, 
and a barrier or turning point in the radial motion is encountered (like in central 
force motion [37]) [38]. The turning point occurs at a non-zero but exponentially 
small radius (for a weak rotational flow) , introducing an additional length- and time- 
scale to the implosion process. Thus, the rotational flow allows us to examine the 
consequences of the singularity formation (and decoupling of the implosion process 
from external scales) on the vibrational mode behavior. 

For the axisymmetric, imploding cylinder, we now allow the liquid to be swirling 
inwards, possessing angular momentum. The velocity field then has an additional 
rotational component which can be written as uq — r/27rr. Kelvin's circulation 
theorem states that the circulation F around a closed loop for potential flow is a 
constant. 



determined by the initial conditions. Thus, the dynamics remains integrable with an 
additional conserved quantity, the circulation. The Hamiltonian in (2.5) must include 
the additional kinetic energy from the rotational flow. This energy is easily found by 
integrating uq over the liquid region and the total kinetic energy is now 




(4.3) 




(4.4) 
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The evolution equation for R{t) (2.4) is modified with an additional term on the right 
hand side, pV'^ /Sn'^R^ . 

The consequence of the rotation is that a singularity in the governing equations 
no longer occurs. Instead, a turning point in the motion occurs when R = which 
can be quickly seen from (4.4) to be exponentially small unless the rotational flow is 
large: Rtum / Roo = exp(— 47ri?/pr^), where we have used the form for the mass from 
Chapter 2, M{R) = 2t^ pR^ hi{Roo / R) ■ This change occurs because the rotational 
flow contains a flnite, and growing amount of energy as the hole collapses, preventing 
the radial velocity from diverging. 

Upon again performing the linear stability analysis for the azimuthal shape per- 
turbations, which now includes terms depend on the rotational flow uq, we flnd a 
description in terms of only the cos(n^) harmonics is insufficient to describe the dy- 
namics. We must now include the second linearly independent sin(n6') harmonics and 
write the perturbed interface as: S{9, t) — R{t) + '^n(^n{t) cos{n9) + Cn{t) sm.{n9)). 
The vibrational modes hn and Cn are coupled by their evolution equations: 

in + (f j 6„ + (1 - „) + „ (^) 6„ + 2n (^) c„ = 0, (4.5) 

c„ + (f ) c„ + (1 - „) + „ (^) J c„ - 2„ (^) 6„ = 0. (4.6) 

For a weak (pF'^/E <^ 1) rotational flow the turning point occurs at an exponen- 
tially small radius, so the leading order radial implosion dynamics are unaffected until 
R — 0{Rturn), outside the limits of a continuum theory. However, the vibrational 
mode behavior does experience a fundamental change - rather than forming standing 
waves with shapes that pulse as i? — > 0, the rotational flow induces traveling waves 
around the circle, again with frequencies proportional to 1/R. Numerical results for 
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bfi and for n = 10 are shown in Fig. 4.3. The fast vibration is the original chirping 
described in Chapter 3 with frequencies y/ri^^/ and the slower oscillation in the 
amplitude envelope is actually caused by the spinning of the quickly pulsing shapes 
as the hole closes. 

To simply see this behavior in (4.5) and (4.6), assume the radial collapse dynamics 
follow R = and uq/uj- = T / na^ <^ 1 . Then trying a priori solutions of the form 

bn = ^0 cos(a;,iln(i?)) cos(-\/n^^^ln(i?)) and Cn = sin(a;,iln(i?)) cos{y/n — lln(i?)) 



yields Un — 




. The full shape can then be written (without the phases) 



S{t,e) = R{t) + ^bn,ocos[uJnHR) - nO]cos{Vn^ln{R)), (4.7) 
n 

clearly showing the traveling wave nature of the solutions, with n-independent fre- 
quencies {uJn/n)R~^ — {V /'Ka^)R~'^ . This is the frequency of the slow oscillations of 
the envelope of the the faster vibrations in Fig. 4.3. Even though the ratio T /na^ 
is not perfectly constant as was assumed in the simple scenario above, but involves 
logarithmic corrections, it changes so slowly that it increased by only a factor of 3, 
over the 10 decades collapse of R. This is reflected in the visible weak variation of 
the long frequency. 

Introduction of the weak rotational flow did not destroy the vibrational mode 
amplitude freezing behavior but did cause the emergence of another frequency in 
the dynamics, the n-independent frequency of the spinning motion of the shapes. 
The surface spins faster and faster as void closes, while the shapes vibrate with 
constant amplitude. The above analysis assumes the rotation is slower than the 
radial converging flow, and thus the effects of the finite turning radius turned out to 
be negligible. It would be interesting to fully examine the asymptotic behavior when 
the rotational flow becomes the same size as the radial flow, and then the destruction 
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Figure 4.3: Coupled mode amplitudes with presence of swirl flow for n = 10, biQ (thick 
line) and ciq (thin line, displaced downwards for clarity). The initial conditions at 
R = Rq are biQ = ciq = bQ with zero derivative. The relative swirl flow strength is 
Uf)/ur ~ 10~^ (over the range of the plot it increases only to 3 • 10~^ as decreases). 
Roo/Rq — 10^ and the implosion energy is about 

of the singularity formation would likely cause rapid amplitude growth or decay. 

In the bubble experiments which start with a quiescent liquid bath, there is likely 
very little circulation present, though it is not possible to completely eliminate. In 
this case, the presence of a weak and unintentional rotational flow are not damaging 
to the vibrational mode behavior, but simply cause a slow spinning of the vibrating 
shapes. Such a rotational flow could be the cause behind observations of rotating 
surfaces and satellite bubbles [39, 40]. 

4.3 Surface tension 



Surface tension forces become negligible to the radial implosion in the asymptotic in- 
ertially dominated regime as seen in (2.5) and verified by experimental measurements 
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of the implosion dynamics. However, here we will see that the effect of surface tension 
on the inertial vibrations within the experimental regime, is to cause a short-lived 
but rapid growth in mode amplitude and increase in frequency in the experimen- 
tal regime before converging onto the asymptotic form given in (3.4). We modify 
(3.4) so that the local curvature is no longer simply 1/R but to first order in hn/R'- 
K pa 4 + (n^ — l)6n cos (n6')/i?^, yielding 



This adjusts the local Laplace pressure jump as we go around the shape through 
the peaks and valleys. When the radial dynamics is inertially dominated, the bn 
coefficient in (3.5) now includes ^n{n -\- l)7rR\a{Roo/ R) / E . By comparison of the 
original inertial term in this coefficient, we see that the cutoff for when surface tension 
will no longer affect the vibrations is when R < R^^^ with the critical radius 



We numerically integrate the solutions for bn{t) using the radial dynamics from a 
numerical solution of (2.7). The solutions show that the effect of surface tension is 
to speed up the vibrations and give additional short-lived amplitude growth when 
the mean radius is large (see Fig. 4.4). Fig. 4.4 (A) shows the growth and frequency 
increase for an n = 10 vibration. The amplitude envelopes for n = 2 — 30 are plotted 
in Fig. 4.4 (B), and the spectrum (using the same definition for the case of viscous 
dissipation above) in (C). Higher modes are impacted more severely and for longer 
into the implosion, since the critical radius Rcrit decreases with increasing n. 

The bubble pinch-off experiments are unique in the fact that they are initially 
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dominated by surface tension (which will efficiently smooth out high modes), but by 
the time the mean radius reaches O(300/im) the radial implosion is inertial. This 
gives surface tension a short time to act on the vibrational modes before the inertial 
limit is reached for R < R^^^ . 
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(a) n = 10 (b) n = 2 - 30 




n 



Figure 4.4: Effect of surface tension on early (large R) vibrational mode behavior. 
Here, Rqo/Rq = 10 and E/tijRq = 10. The amplitudes are rescaled by their initial 
values 6o the mean radius is rescaled by its initial value Rq. (a) Representative 
vibration for n = 10 (thin line) shows that surface tension causes an initial increase in 
the amplitude (envelope is thick line) and frequency of vibration. After the mean ra- 
dius has decreased beyond R^q^ ~ Rq/IO'^, the vibration settles in to the asymptotic 
inertial state with amplitudes and periods as in Fig. 3.1. (b) Envelope of vibrational 
mode amplitudes for n = 2 — 30 (lower- to upper- most), (c) Amplitude spectrum 
taken at R/Rq — 10~^ (open circles) and at R/Rq — 10~^ (solid circles). 



CHAPTER 5 
CONNECTION BETWEEN THEORY AND 
EXPERIMENTS 

The preceding chapters laid out a scenario of how and why a nearly cylindrical void 
collapsing in an inviscid liquid retains an imprint of initial distortions to its shape. 
The cylindrically-symmetric model of Chapter 2 predicted how the void radius should 
collapse to a point, and the stability analysis of Chapter 3 predicted that and initial 
distortion causes vibrations of the shape as the implosion proceeds. The vibration 
amplitudes of all wavelength disturbances become frozen, and eventually overtake the 
leading-order symmetric implosion, generating strong distortions. Consideration of 
viscous dissipation, rotational flow, and surface tension effects in Chapter 4 predicts 
that the mode amplitude freezing and chirping (frequencies oc 1/r) are characteristics 
within the realm of experimental observation. 

Below we briefly review two types of experiments that generate void collapse in 
water. Then, because the vibrational behavior is inherently linked to the radial im- 
plosion dynamics via (3.4) we show the quasi-2D is capable of describing the evolution 
of 3D experimental profiles. Finally, we directly compare an experimentally observed 
shape vibration to the prediction of (3.4) for an n = 2 vibrational mode. 



5.1 Review of experiments 

In the bubble release experiments a bubble is grown quasi-statically by slowly pump- 
ing gas through a nozzle [9, 10, 11]. Surface tension keeps the bubbles connected 
to the nozzle until it can no longer stabilize the surface against the upwards pull of 
buoyancy. At this point, the neck of the bubble rapidly pinches off. The entire dis- 
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connection process occurs on a ms time scale, much faster than the rising motion of 
the bubble. The earliest dynamics are determined by a competition between surface 
tension and gravity. However, once the implosion begins, the inertia of the water 
rushing in to fill the void quickly starts to dominate the process, evidenced by the 
time evolution of the minimum radius /i^riin following close to a 1/2 power law in 
time, rather than 2/3 if surface tension were to drive the flow. 

For the disc impact experiments, a solid disc is rapidly drawn into a bath of 
water, creating a cavity in its path [14, 12, 16]. Because of the high impact velocities 
used, the initial Weber number, which characterizes inertial effects relative to surface 
tension is greater than 10^ {We = pU'^h/'j] p is the liquid density, 7 is the surface 
tension of the gas/liquid interface, U a characteristic velocity scale, h a characteristic 
length scale). Thus there is no regime here where surface tension can efficiently act 
to smooth out distortions to the cavity shape. 

Both types of experiments observe a medley of different kinds of disconnection 
events depending on the specifics of the set-up, geometry and orientation of the nozzle 
or disc, and control parameters. In the bubble experiments, side-by-side satellite 
bubbles remain after disconnection if the nozzle is not carefully leveled [9, 10]. Fig. 
1.1 (C) shows the bubble neck before and after disconnection when the nozzle has been 
purposely tilted by 2°. Even more exotic and highly asymmetric shapes are observed if 
the bubble is not grown quasi-statically from the circular nozzle, but is instead blown 
rapidly through a slot shaped nozzle [10]. Often times in the impact experiments, 
the surface before disconnection appears textured, and after disconnection leaves 
behind a disordered swarm of satellite bubbles, rather than cleanly disconnecting at 
a point [14, 12]. These observations signal an instability in the underlying dynamics 
to azimuthal and/or axial perturbations, which prevents the axisymmetric implosion 
from being realized under generic conditions. Similar effects have not been observed in 
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experiments where surface tension is the force driving breakup. For example, a water 
drop breaking in air appears to always break in a universal way, independent of the 
exact experimental set-up [18]. Despite the clear evidence for the initial distortions 
in the azimuthal direction playing a critical role in the dynamics before and after 
bubble disconnection, it remained unexplored in the context of cylindrical voids until 
recently [13]. The accumulation of experimental results suggests that the fixed point 
for axisymmetric cavity collapse proposed [8] is not stable once the azimuthal direction 
is considered. 

The final stages of both implosion situations have the same inertia- dominated dy- 
namics, but the initial states are quite different. The impact experiments never pass 
through a surface tension dominated regime, whereas the bubble experiments always 
start there. This puts both experiments in unique positions. For impact, large am- 
plitude and/or short wavelength distortions can be applied by using different shaped 
discs, because surface tension will not smooth them out. However, the advantage in 
the bubble experiments is the possibility to generate smooth, low wavelength, single 
mode perturbations to the shape. 

5.2 Axisymmetric implosions 

With the implosion dynamics outlined above we can now compare the profiles pre- 
dicted by conservation of energy to those actually observed in experiments. Neck 
profiles from traces of experimental images of a bubble released from a 4 mm diame- 
ter circular nozzle are shown in Fig. 5.1 (b). To generate the theoretical profiles for 
comparison we first use the time-dependence of the minimum (at 2; = 0) to define the 
implosion energy per unit length. Within the range shown, the radius is following 
close to a 1/2 power law, signaling that the radial implosion is now far within the 
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Figure 5.1: (a) Collapse of minimum radius for air bubble released from circular, 
4mm diameter nozzle (red circles). Theoretical inertial implosion from (2.7) with an 
implosion energy per unit length of = 12 dyne (black line). Dashed line shows 
slope of 1/2. (b) Experimental neck profiles of the same air bubble (blue diamonds). 
Theoretical profiles calculated as described in text using conservation of energy for 
each isolated layer (lines). Profiles are at 4.5, 49, 108, 166, 270, and 372 /is before 
disconnection (inner to outer curves). Experimental data provided by Nathan Keim 
(U. of Chicago). 
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inertially dominated regime. We then compare the evolution of the minimum to that 
predicted in the inertial regime, from numerical integration of (2.7) to determine the 
implosion energy E. 

Taking Rqo = 2mm, a physically relevant length scale as it is the radius of the 
nozzle, we find E — 12 dyne gives the best fit to the minimum (Fig. 5.1 (a)). The 
result is insensitive to the value of Rqo used and good fits within experimental error 
can be found for R^q = 1mm with E = 8.7 dyne or up to Rqq = 4mm with E = 15 
dyne. These results would be indistinguishable by eye compared to those shown in 
Fig. 5.1(a). 

The implosion energy of 12 dyne was then applied at every height to track the 
full profile evolution. Because of the time invariance in the problem, the symmetric 
form of the profile closest to disconnection (innermost profile, averaged about z = 0) 
was used as the initial state, and then the expansion of the radius at each height 
was evolved according to R — +•>/ E/i^pR^ \n{Roo/R)- Though the angles observed 
in bubble experiments typically are not small (half angles range from about 25-35 
degrees depending on conditions), the quasi-2D model does a surprisingly good job 
describing the complete 3D surface evolution. This suggests the errors involved in 
applying the quasi-2D model are quite small over the finite range in /imm afforded 
by the experiment, though this remains to be checked via numerical simulation. 

One aspect of the experiments which is quite obviously not captured in the model 
is the development of an up/down asymmetry in the cones above and below the 
minimum. The upper cone has an opening angle which is larger than the lower 
by about 5 degrees, but the precise amount depends on the vertical range used to 
define the angle. One might think that this could be simply accounted for within 
the framework of the 2D model by incorporation of a hydrostatic pressure gradient, 
which would break the axial symmetry about the minimum. 



41 

To introduce the hydrostatic pressure gradient —Apg we simply include it in the 
Hamiltonian (2.5) so the pressure term becomes (AP — Apgz)7rR'^. Then, the energy 
is a function of z, with deeper layers having larger potential energies because the 
hydrostatic pressure at r = Rqq is higher there. If the radius at the same height 
above and below z — are equal, the one with higher energy has a faster implosion 
rate (2.7). Thus, naively, one might expect that the cone below the origin will implode 
faster and have a smaller angle relative to the upper cone. However, one finds that 
this is not the case. Instead of causing a change in the cone angles, the primary effect 
is a downwards linear translation of the profile with the angles unaffected. 

The profiles can be shifted to follow this translation of the minimum and are of 
the form R(ijj)'^ — h"^^^ + 2cqHq'iIP' + C2'0^ where h^i^ is the radius of the minimum 
at 2; = ^TOin s-iid = z — Zjj^in- The cone angles remain \/2cqHq. The only possibility 
for the observed difference in apparent cone angles is if the if)^ term is very large, 
which requires the hydrostatic energy to be a substantial component of the implosion 
energy. This case is relevant to the larger-scale impact experiments [12, 14] but not for 
bubble pinch-off where the scale of the neck is O(100/xm). An additional inconsistency 
is fact that if the difference in cone angles was caused by a large hydrostatic pressure 
gradient, the motion of the minimum would be downwards, but the experiments 
observe and the upwards motion of the minimum. Within the framework above, the 
hydrostatic pressure alone is not sufficient to explain the difference in cone angles 
and upwards motion of the minimum observed in the experiments. The up/down 
asymmetry which appears may be partly due to the slow upwards motion of the 
buoyant bubble, but is likely related to non-local flow coupling or energy transfer in 
the vertical direction, and is a point to be examined further in the future. 
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5.3 Predicted and observed vibrational mode 

Recent experiments at the University of Chicago have been able to observe vibrations 
of the bubble neck shape as it implodes by releasing a bubble from an oblong noz- 
zle [13, 39]. As the bubble grows from the nozzle, the neck is initially elongated in 
the same direction as the slot. Though surface tension is efficient in rounding up the 
bubble neck, it does not completely erase the effect of the highly asymmetric nozzle, 
and a small asymmetry in the cross section can be seen. 

As the neck begins to rapidly implode, the neck shape vibrates so that at a later 
time it will be elongated in a direction opposite to slot. Depending on the speed 
at which the bubble is released from the slot nozzle (either quasi-statically or via a 
rapid burst of air) different numbers of periods can be observed. Fig. 5.2 shows these 
vibrations for quasi-static release from a slot of size 1.6 mm by 9.6 mm. By recording 
the implosion with two high speed cameras at a 90 degree angle, the perturbation 
amplitude is simply the half the difference in the measured widths from the two views 

02 = (hside - h front) 1'^^ and the mean radius is ^ = {hgide + hfront)/'^- 

Above a measured mean radius of ^ ~ 400//m there is a rapid decay of the 
distortion amplitude which is likely due to the fact that the dynamics is truly three- 
dimensional and dominated by surface tension in this regime. Around R ~ 300//m 
a vibration centered around a2 = sets in and continues until the last possible data 
point when R ^ 20/xm. Below R ~ lO/^m there is evidence of some weak decay in the 
amplitude, which could be due to viscous dissipation. Though the vibration shown 
in Fig. 5.2 is at height of the minimum radius in the experiment, it is also possible to 
observe the vibrations vibrations at different heights along z. 

Below we outline the parameters, which are consistent with the experiment, used 
in the model to generate the theoretical vibration shown in Fig. 5.2. We begin by 
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Figure 5.2: Experimentally observed vibrations of n = 2 mode by release from 1.6 
mm by 9.6 mm slot nozzle (red). Theoretical vibration calculated from (4.8) as 
described in text (solid black line). Inset: Experimental mean radius collapse (red) 
corresponds to an inertial implosion with £" = 12 dyne, and Rqq = 5mm (solid black). 
A power-law of 1/2 is shown by the dashed line for comparison. 

defining the total energy of the symmetric radial collapse which reproduces the radial 
inertial implosion dynamics for R{t). For the data shown in the inset to Fig. 5.2, this 
value is = 12 dyne/cm when i?oo = 5mm, the average radius of the slot. Similarly 
to the data in Fig. 5.1, while Rqq is chosen within a physically meaningful range 
as the system size, the radial dynamics from energy conservation (2.7) describe the 
observed trajectory. Below we estimate the appropriate value for 7 to be used in the 
evolution equation for a2{t) (4.8) which will take into account the axial curvature 
contribution to the Laplace pressure jump across the interface as done in (2.16). 

From examination of experimental images of the data set used to generate the 
data shown in Fig. 5.2 we find that -Raa;/^mMi ~ 2.8. The value of 2.8 was chosen 
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by taking the average of the ratio Rax/hmin over the range of the experimentally 
observed vibration. Thus, to effectively include the axial contribution, a value of 
7 = (1 — 1/2.8)70 = 0.6470 = 46 dyne/cm where 70 is the surface tension of an air 
water interface at room temperature, 72 dyne/cm, was used in the numerical solution 
of (4.8). The black hne plotted in Fig. 5.2 is the result of a numerical solution to 
(4.8) with the parameters outlined above. The amplitude and phase (not predictable 
in this model) were chosen to fit the peak around R = 100/xm. 



CHAPTER 6 
CONCLUSION 

The 2D model of the implosion of a cylindrical void in water is an integrable Hamil- 
tonian system with special properties. These properties include conservation of the 
implosion energy and an exceptional perturbation mode spectrum which freezes as 
the implosion proceeds, and causes constant amphtude inertial vibrations about the 
radial implosion. This ideal system is closely related to the dynamics of the thin neck 
of air which collapses as an air bubble disconnects from an underwater nozzle. The 
memory mechanism in the simplified model is capable of explaining why asymmetries 
are observed to persist in experiments, and has predictive power for how azimuthal 
disturbances on the surface can dramatically alter the final stage of the pinch-off 
process. 

The main features (amplitude freezing and vibrational chirping) of the inertial 
vibrations are robust against additional effects of surface tension, viscous dissipation, 
and rotational flow within the experimental regime, and thus we are able to directly 
compare the predicted vibrational mode to one excited in experiments. However, 
when comparing the observed amplitude and frequency of the n — 2 vibration, the 
most important effect was found to be surface tension. Inclusion of surface tension 
effects via the Laplace pressure jump across the interface acted to speed up the 
vibrations for a short time within the experimental regime. 

At this point, a variety of physical effects were considered separately, to isolate 
the individual consequences on the vibrational mode spectrum. There is the potential 
for interesting mode spectrums to arise if multiple effects are combined. For example, 
the form of the normal stress balance on the interface will change substantially if 
both viscous dissipation and a rotational flow are included. It may also be possible 
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to find peaks in the spectrum, where one wavelength grows more rapidly than others, 
when combining surface tension (which causes amplitude growth) with viscous dis- 
sipation (which causes amplitude decay). More exploration of the outcome of these 
possibilities will be necessary to predict vibrational mode behavior for a larger range 
of experimental parameters. 
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